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Abstract 

The quantum coupling of fully different degrees of freedom is a challenging path towards new 
functionalities for quantum electronics m- Here we show that the localized classical spin of a 
magnetic atom immersed in a superconductor with a two-dimensional electronic band structure 
gives rise to a long range coherent magnetic quantum state. We experimentally evidence coherent 
bound states with spatially oscillating particle-hole asymmetry extending tens of nanometers from 
individual iron atoms embedded in a 2H-NbSe2 crystal. We theoretically elucidate how reduced 
dimensionality enhances the spatial extent of these bound states and describe their energy and 
spatial structure. These spatially extended magnetic states could be used as building blocks for 
coupling coherently distant magnetic atoms in new topological superconducting phases mn]. 
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Coupling different degrees of freedom of both quantum and classical objects yields new 
quantum functionalities not available in each system taken separately. In this regard, new 
hybrid quantum systems have been recently designed such as individual atoms and optical 
cavities coupled through photon exchange [Tj , or a single quantum dot coupled to a mechani¬ 
cal oscillator via strain [2j. A remarkable example of emerging phenomena is the observation 
of Majorana end states resulting from the coupling of a ferromagnetic chain of Fe atoms and 
a superconducting substrate |3j. These states are in the focus of numerous theoretical works 
in the held of quantum computing as they are associated with non-abelian statistics. Yet, 
the spatial extent of Majorana end states measured by S. Nadj-Perge et ah jTj is restricted 
to a few atomic distances, making difficult to handle them for braiding. 

An alternative proposal for manipulating Majorana quasiparticles consists in engineering 
a one-dimensional topological superconductor in a chain of magnetic atoms with a spiral 
magnetic order on the surface of a superconductor HHII]. Individual local magnetic mo¬ 
ments act destructively on Cooper pairs, leading to discrete spin-polarized states inside 
the superconducting energy gap, predicted by Yu, Shiba and Rusinov [T3IIT5] (YSR). Rusi- 
nov suggested that around magnetic atoms the decaying YSR wavefunction should have a 
spatially oscillating structure naiisi [16]. The emergence of topological superconductiv¬ 
ity depends on the YSR states mediated coupling inside the magnetic chain m- While 
Nadj-Perge et al. j3] relied on a direct interaction between neighboring atoms to generate 
short-ranged Majorana quasi-particles, in the latter case the characteristic length is that of 
the YSR bound states, which may extend up to the scale of the superconducting coherence 
length. Enhancing the spatial extent of YSR bound states would facilitate the remote cou¬ 
pling of magnetic systems through a superconducting state, opening the route towards an 
easier manipulation of Majorana quasiparticles and the creation of new topological quantum 
devices. 

Here we reveal that the dimensionality plays a critical role in the spatial decay of the YSR 
bound states. Calculations using Rusinov’s approach na are presented in Figs. [^A-D. They 
show that a three-dimensional isotropic s-wave superconductor induces a dramatic decay of 
the YSR states out of the magnetic impurity. These calculations are in agreement with the 
atomically short spatial extent of the YSR states observed in all previous scanning tunneling 
microscopy (STM) studies of single magnetic impurities in superconductors: Co, Cr, Mn, 
Gd atoms deposited on Pb or Nb crystals dZlIIB] and manganese-phtalocyanine molecules 
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FIG. 1: Comparison between 3D and 2D for the spatial extent of Yu-Shiba-Rusinov 
states (A) to (D) ((E) to (H)) Calculated behavior of a Yu-Shiba-Rusinov bound state in an 
isotropic s-wave superconductor with three-dimensional (two-dimensional) electronic band struc¬ 
ture. (A) and (E) are schematic views of the interaction of Cooper pairs with a classical magnetic 
impurity. (B) and (F) are calculated scanning tunneling spectra at various distances from the im¬ 
purity showing the fully polarized YSR states inside the superconducting gap. (C) and (G) are 
simulated conductance maps around the impurity showing the spatial extent of one peak of the YSR 
state presented in Figs. B. and F. respectively. (D) and (H) are simulated conductance between 
-0.6 and 0.6 mV along the dotted line out of the impurity in Figs. (C) and (G) respectively. 

deposited on Pb [12]. By extending this theory to two dimensions (2D), we evidence in Figs. 
[^E-G that superconductors with two-dimensional electronic structure should host YSR 
bound states with spatial extent orders of magnitude larger. Hence layered superconducting 
materials such as 2H-NbSe2, known for their two-dimensional character, are good candidates 
for supporting these long range quantum states |20j . 

In this work, we studied single crystals of 2H-NbSe2 containing a few tens of ppm of mag¬ 
netic Fe impurities. Atomically resolved topographic STM imaging of the surface (Fig. j^B) 
exhibits the characteristic charge density wave pattern. Individual impurities appear as 
bright spots. Whereas on topographic STM images the magnetic and non-magnetic im¬ 
purities cannot be distinguished, the magnetic impurities are the only ones to present a 
characteristic spectroscopic signature inside the superconducting gap (see supplementary 
SI). 

The scanning tunneling spectroscopy studies performed at 320 mK, well below the critical 
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FIG. 2: Structural and superconducting properties of 2H-NbSe2 (A) Atomic structure of 
a 2H-NbSe2 crystal. (B) 19x17 nm^ topographic image of a 2H-NbSe2 sample with the atomic 
lattice modulated by a charge density wave. The image is taken at V=-200 meV and 1=80 pA. (C) 
Abrikosov lattice in 2H-NbSe2 showing the star shaped structure of vortices in a magnetic held of 
0.1 T. (D) Conductance map taken at V=-0.05 meV showing a few star-shaped structures created 
by localized magnetic impurities at zero magnetic held. Measurements were performed at 320 mK. 

temperature of 7.2 K, reveal YSR bound states around the randomly dispersed magnetic iron 
impurities in 2H-NbSe2 (see supplementary S2). Remarkably, these states are characterized 
by a six-pointed star shaped electronic signature extending as far as 10 nm from defects, 
as can be seen in Fig. iD. This is more than 10 times larger than the previously observed 
extension of a few A for YSR bound states [IM9] and is comparable to the in-plane coher¬ 
ence length of 2H-NbSe2. This long range pattern is due to the two-dimensional character of 
2H-NbSe2 (see Fig. [^, compared to the observed short range in three-dimensional materials 
such as Pb or Nb. This unusually long spatial effect may also be amplihed by the fact 
that the Fe impurities are embedded in the atomic lattice. Therefore they may experience a 
stronger electronic coupling to the superconducting condensate than the adsorbed impurities 
used in previous experiments. 

The arms of the star shaped pattern of YSR states are turned by 30° with respect to the 
crystallographic axes of 2H-NbSe2 (Fig. [^A), which corresponds to the reciprocal lattice 
vectors (a* and b*). This orientation is also the same as the one of the star-shaped vortices 
observed in 2H-NbSe2 (Fig. [^C). In addition to the dominant type of impurities shown on 
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Fig. §D and [^A we also observe impurities of Cr and Mn giving star shaped structure 
with the same orientation but with a thicker pattern and slightly different YSR energies 
(see supplementary S4). Therefore this six-fold symmetry is likely to arise from a common 
origin, and in both cases reflects the anisotropy of the Fermi surface jlO], as supported by 
our simulations. 

The tunneling spectra acquired over a chosen Fe impurity (see spectroscopic map in 
Fig. [^A) show a YSR bound state which takes the form of two peaks at positive and 
negative energies {Esuba — ± 0.2 A) inside the superconducting gap of 2H-NbSe2 (red 
curve in Fig. |^B). Apart from the YSR state the characteristic superconducting spectrum 
is perfectly preserved. The YSR peak at negative bias is much stronger than the one at 
positive bias, highlighting a strong particle-hole asymmetry near the magnetic atom, as 
presented in Figj^H. The presence of a single pair of YSR peaks in the gap indicates that 
s-wave diffusion channel {I = 0) dominates, which suggests that the iron impurities may 
be considered as punctual defects. Similarly, in Gd/Nb, Mn/Nb im or Mn-Phtalocyanine 
on Pb IIS] the I = 0 diffusion channel was the only one to be activated. On the contrary 
for Mn/Pb, YSR states for I = 0 and I = 1 were observed and I = 2 states were found 
for Cr/Pb(lll) [18]. However in all these works, the spectroscopic signatures associated 
with the impurities completely vanished a few A from their center. In this context our 
measurements show that the local nature of the interaction does not prevent the existence 
of long range effect on the density of states. 

In order to recover the symmetry of the observed YSR bound state one needs to take into 
account the band structure of the material. The hexagonal symmetry observed experimen¬ 
tally in Figs. iD andj^A is well reproduced in the framework of the Bogoliubov-de Gennes 
formalism [21]. This is done by numerically solving the Schrodinger equation with an almost 
exact tight-binding description of the band structure of 2H-NbSe2 (see supplementary S5). 
As we only observe I = 0 states we assume a strictly on-site interaction while treating the 
magnetic impurity classically, i.e. assuming a large spin number S (see supplementary S3). 
The interaction potential contains both a magnetic and non-magnetic part and reads as 


JS 


Hjmp — „ ^04, A A"(cQ.j.Co'|- -|- Cgj^Coj,) 


( 1 ) 


Where the cq and cj, operators are respectively the annihilation and creation operators for 
electrons with spin a on the magnetic atom site. The hrst term corresponds to a Zeeman 
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splitting between spin up and spin down electrons for a coupling strength J/2 between the 
superconducting electrons and the individual atom. The second term is a non-magnetic 
diffusion potential of amplitude K. Using this approach, we recover the typical star shaped 
structure as presented on Fig. [^aligned along the reciprocal lattice vectors. 




FIG. 3: Spectral and spatial properties of an extended Yu-Shiba-Rusinov bound state 
in 2H-NbSe2 (A) Experimental conductance map taken at -0.13 meV. Two a and b lines indicate 
the crystallographic axes of 2H-NbSe2 while the a* and b* indicate the directions in the reciprocal 
space. (B) Characteristic experimental spectra taken on top of the impurity (red), on the left 
branch, 4 nm from the center of the impurity (green) and far from the impurity (blue). (C) Spatial 
and energy evolution of the experimental tunneling conductance spectra, dl/dV{x,V) along one 
branch of the star. The left side of the figure corresponds to the center of the star and the right 
side to the top-right corner of the scanning area. The color scale is the same as the one used in 
(B). (D) Conductance profile of the electron and hole like YSR states as a function of the distance 
to the impurity along the same line as for (D). 


We now focus on the details of the star presented on the experimental spectroscopic map, 
taken at the energy of the strongest YSR peak -0.13 mV (see Fig. [^A). The high tunneling 
conductance in the center of the star (red color) corresponds to a very strong peak in the 
tunneling spectra (red curve on Fig. [^B) localized on the impurity. On the surrounding 
tail the amplitude of YSR peaks decreases as shown by the green conductance curve of 
Fig. [^B acquired at 4 nm away from the impurity. This decrease is oscillatory, resulting in 
interference fringes with a periodicity of 8 A clearly visible on the conductance map. The 
evolution of the conductance spectra along one arm of the star is shown in Fig. [^C. The 
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interference fringes in the conductance for the electron and hole excitations are in an almost 
perfect spatial antiphase. It appears that a few nanometers away from the impurities the tail 
of the YSR states exhibits a similar amplitude for the electron-like and hole-like excitations 
(see Figs. gC^^D), i.e. the particle-hole symmetry is progressively restored away from the 
defect. Our analysis shows that the oscillations are due to the scattering of electrons by the 
impurity at the saddle points between the pocket around F and the pockets around the K 
points of the Fermi surface of 2H-NbSe2. 

The observed oscillations being of the order of the Fermi wave-length, they cannot be cap¬ 
tured by a discrete tight binding model. Following Rusinov m and taking the continuum 
limit description of the equation giving access to the YSR states, we extract the asymp¬ 
totic behavior of the wave function far from the magnetic impurity (see supplementary S7). 
We qualitatively reproduce the two characteristic length scales of the experimentally ob¬ 
served interference pattern. Assuming an isotropic energy band, the YSR energy can be 
parametrized as E = Acos((5’'' — 5“) jT^] with tan 5^ = Kuq ± JS/2vq where z/q is the den¬ 
sity of states at the Fermi energy. In two dimensions the YSR wave function can be shown 
to behave at large distance from the impurity as (see supplementary S7): 

V-±(r) = j (2) 

where 4^+ V’- denote respectively the electron and hole components of the YSR wave 

function -0, N is a normalization factor, is the Fermi wave vector and vp the Fermi 
velocity. This behavior presented in Fig. [^H is in excellent agreement with the experimental 
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FIG. 4: Tight binding calculation of the spatial structure of Yu-Shiba-Rusinov states 

Local density of state (LDOS) computed with JS/2 = 120 meV and K = —180 meV. a) LDOS for 
electron-like YSR state with energy Eghiba = 0.59A. b) LDOS for hole-like YSR state with energy 
Eshiba = -0.59 A. 
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Figs. [^C-|^D. This result highlights the dimensionality dependence as the decay of the local 
density of states goes as 1/r in 2D and 1/r^ in 3D (see supplementary S8). Furthermore, 
for deep YSR states, which corresponds to a dephasing — 6~ —> ±7r/2, the electron 
and hole YSR states are indeed in antiphase far enough from the impurity. In our approach 
each component of the YSR states is fully polarized both in spin and charge. The analytical 
equation links the Friedel-like long distance decay of the YSR state to the superconducting 
coherence length and the oscillatory behavior to the Fermi wavelength. 

In conclusion, by coupling a classical spin to a superconductor with a two-dimensional 
electronic structure we were able to unveil a long range coherent magnetic quantum state 
with a spatially oscillating electron-hole asymmetry. As this effect is related to the dimen¬ 
sionality it should manifest in a wide variety of superconductors, such as lamellar materials 
or recently discovered superconducting monolayers of Pb/Si(lll), In/Si(lll) |2^ |23] and 
FeSe/SrTiOaIJS]. The interaction between long range YSR states has now to be explored. 
Then, it could be used for producing new topological phases in hybrid systems. Arrays of 
magnetic atoms and molecules coupled trough a superconducting medium are indeed ex¬ 
pected to present a large variety of topological orders. For instance, a chain of magnetic 
atoms coupled through the spatially extended YSR bound states with a helical spin order 
could lead to a topological triplet superconductivity with Majorana quasiparticles at its 
extremities IMI]. Such spiral magnetic ordering could be stabilized using the RKKY inter¬ 
action which oscillates with inter-atomic distance, as for Co/Ir(001) [ 22 ]- The fine tuning of 
the YSR and RKKY mediated interactions could be achieved using one- or two-dimensional 
self-assembly of atoms and molecules on templated substrates P EH EHl ET]. 
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